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We have searched for CP violation in the charmed meson decays D+ s) -> K%n+ and D+ a) -> K° S K+ 
using 673 fb -1 of data collected with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. No evidence for CP violation is observed. We report the most sensitive CP asymmetry 

measurements to date for these decays: A^, p ^ KsTr — (—0.71 ± 0.19 ± 0.20)%, A p^ s * = 

(+5.45 ±2.50 ±0.33)%, A^^ K ° sK+ = (-0.16 ± 0.58 ± 0.25)%, and = (+0.12 ± 0.36 ± 

0.22)%, where the first uncertainties are statistical and the second are systematic. 
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Violation of the combined Charge-conjugation and 
Parity symmetries (CP) in the standard model (SM) 
is produced by a non-vanishing phase in the Cabibbo- 
Kobayashi-Maskawa flavor- mixing matrix For 
charged meson decays this may be observed as a non- 
zero CP asymmetry, defined as 



A C p = 



T{X+ -+ /+) - T(X- -+ /") 

rpr+->/+) + rpr-->/-) 



(i) 



where T is the partial decay width, X denotes a charged 
meson, and / is a final state. 

In the SM, the charmed particle processes for which 
a significant non- vanishing CP violation is expected are 
singly Cabibbo-suppressed (SCS) decays in which there is 
both interference between two different decay amplitudes 
and a strong phase shift from final state interactions. In 
the SM, CP violation in SCS charmed meson decays is 
predicted to occur at the level of O(0. 1)% or lower 
The SM also predicts a CP asymmetry in the final states 
containing a neutral kaon that is produced via K° — K° 
mixing even if no CP violating phase exists in the charm 
decay amplitudes. The expected magnitude for this type 
of asymmetry is (0.332 ± 0.006)% from if£ semileptonic 
decay Q . Searches for C P violation in charmed mesons 
are complementary to those in B and K mesons, since the 
former tests the CP violating couplings of the up-type 
quarks while the latter those of the down- type quarks. 

In this Letter we report results from searches for CP 
violation in the D+ ->■ K s ir+, D+ ->• K° S K+, D+ ->• 
K s n + , and D+ -> K S K + decay processes [3]. The for- 
mer two channels are SCS decays and the latter two are 
mixtures of Cabibbo-favored (CF) and doubly Cabibbo- 
suppressed (DCS) decays, where SM CP violations de- 



scribed above are expected. Observing Aqp values of 
0(1)% or larger in the decays considered in this Letter 
would represent strong evidence for processes involving 
physics beyond the SM 

The data used in this analysis were recorded at or near 
the T(45) resonance with the Belle detector Q at the 
e + e~ asymmetric-energy collider KEKB The sample 
corresponds to an integrated luminosity of 673 fb -1 . 

D + and Df candidates are reconstructed using the 
same requirements used in the measurement of the 
branching ratios for these same decays reported in 
Ref. Q. Figure [T] shows the reconstructed K^tt^ and 
KgK^ invariant mass distributions. All signals are pa- 
rameterized as a sum of two Gaussian distributions ex- 
cept for Df — > KgTT + in which a single Gaussian is used 
for the signal parameterization. The parameterizations 
of the random combinatorial background and the peaking 
background due to K/ir misidentification are described 
in detail in Ref. [8|. 

X + ^-K°h + 

We determine the quantities A cp s defined in 
Eq. ([1} by measuring the signal yield asymmetry 



A 



>K°h+ 



N x 

Jv rcc 



X + ^K%h+ 



n; 



x- 



(2) 



where AT reG is the number of reconstructed decays and 
ft is a charged hadron. The measured asymmetry in 
Eq. ^ includes two contributions other than Aqp- One 
is the forward-backward asymmetry (Afb) due to 7* — Z° 
interference in e + e~ — > cc and the other is a detec- 
tion efficiency asymmetry between positively and nega- 
tively charged tracks (A^ ) . Since Kg mesons are recon- 
structed from a tt + it~ pair, there is no detection asym- 
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each bin. Figure[2]shows the Acp map of D + — > KgTr + in 
bins of (p l £ h , cos0 l £ h , cos6>^+ s ). Calculating a weighted 

average of the Acp values over the 3-D bins, we obtain 

Ac + p^ K ° s7T+ = (-0.71 ± 0.26)% where the uncertainty 
originates from the finite size of the D + — > KgTr + (0.19%) 
and D+ -> c/>7r + (0.18%) samples. The x 2 /d.o.f over the 
3-D bins is found to be 31.4/24 which corresponds to 14% 
probability. 

The statistical precision of the — > KgTT + sample 

is too low to allow for a 3-D correction to A IC c s 
For this mode we correct for asymmetries other than 
Acp with an inclusive correction obtained by subtract- 



ing Arec~"^ s " from A^ p ^^ s " after integrating over 
the entire (p^ b , cos#Jf b , cos9™ s ) space. The inclusive 
correction is (— 0.34±0.18)% where the uncertainty is en- 
tirely due to the statistical uncertainty of the D+ — > <frir + 

sample. The value of A CP s is measured to be 
(+5.45±2.50)%, where the uncertainty is statistical only. 



FIG. 1: Invariant mass distributions for the Kgn ± and KgK ± 
final states. Points with error bars (note the small size of them 
due to the large sample) show the data and the histograms 
show the results of the parameterizations of the data. Signal, 
peaking background, and random combinatorial background 
components are also shown. 
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metry other than A^ 
as 



Eq. can therefore be expressed 



A 



A 



CP 



A x+ 

A FB 



A 1 ; 



(3) 



To correct for the asymmetries other than Acp, we use 
reconstructed samples of Df — > (f>n + and D° — > K~n + 
decays and assume that Aqp in CF decays is negligibly 
small at the current experimental sensitivity and that 
Apb is the same for all charmed mesons. We reconstruct 
(f> mesons via their K + K~ decay channel for Df — > c/>7r + , 
requiring the K + K~ invariant mass to be between 1.01 
and 1.03 GeV/c 2 . 

The measured asymmetry for Df — > (f>Tt + is the sum 

D + + 

of A F S B and A^ . Hence one can extract the Ac p value 
for the Kg7r + final states by subtracting the measured 
asymmetry for Df —> <pTT + from that for — > KgW + . 
The subtraction is performed in bins of tt + momentum, 
p^ b , and polar angle in the laboratory system, cos6*Jf b 
(because A^ + depends on these two variables while it is 
uniform in azimuthal angle), and the charmed meson's 
polar angle in the center-of-mass system, cos #^ms ( smce 



cos9™° is correlated with cos#if b and A F { g depends on 



it) . The choice of the three-dimensional (3-D) binning is 
selected in order to avoid large statistical fluctuations in 
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FIG. 2: Measured Acp values for D + 

(p^ h , cos6^ h , cos6'^+ s ). Empty bins where no entries are 
plotted have no statistics. 



The dominant source of systematic uncertainty in 
the A C p measurement is the uncertainty in the 

D~^~ — xh-K^ 

And determination, which originates from the fol- 

lowing sources: the statistical uncertainty of the selected 
D+ (fm + sample (0.18%); the choice of the M(K + K~) 
interval (0.03%); and the choice of binning for the 3- 

D map of A-ec (0.03%). Another source is the 

choice of fitting parameters for the invariant mass dis- 
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tributions: binnings, mass windows, and background 
parameterizations together, contribute uncertainties of 



0.04% to Acp^^ s7r ' and 0.27% to A^" s " , where 

D + ->K°-k + 

the larger uncertainty on A CP s is inherited from 
the low statistics of Df —> K^tt + . We also consider pos- 
sible effects due to the differences in interactions of K° 
and K° mesons with the material of the detector. K a 
and K° mesons considered in this Letter are produced 
through the weak interaction and interact with the ma- 
terial near the interaction point until they decay into 
7r + 7r~ pairs. This produces a non-vanishing asymmetry 
originating from the different strong interactions of K 
and K° mesons with nucleons. Assuming that the dif- 
ferences between K° and K° interactions with nucleons 
are the same as those for K + and K~ interactions, we 
calculate the probability of K° and ^"-nucleons interac- 
tions using the known K + and K~ cross sections Q and 
take into account the time evolution of neutral kaons. 
We consider the beam pipe [6] and the silicon vertex de- 
tector in our estimates of the i^/J^-material effects. 
The uncertainty in the CP asymmetry due to K°/K°- 
material effects is found to be 0.06%. A summary of sys- 



tematic uncertainties in A C p 



is given in Table fl] 
By combining all systematic uncertainties in quadrature, 



we obtain A 

A CP 



CP 



(-0.71 ± 0.19 ± 0.20)% and 

(+5.45 ± 2.50 ± 0.33)%, where the first 
uncertainties are statistical and the second are system- 
atic. 

The method for the measurement of Aqp in the K^K + 
final states is different from that for the Kgir + final 

states. The A F £ and A™ components in A^c 

are directly obtained from the Dj — > (f>Tr + sample, but 

there is no corresponding large statistics decay mode that 

can be used to directly measure the A F £ and Af com- 

ponents in A rcc B • Thus, to correct the recon- 

structed asymmetry in the K { gK + final states, we use 
samples of D° —> K~ir + as well as — > <fiir + decays. 

The measured asymmetry for D° — > K~tt + is a sum 
of Ap° B , Af, and A* + . Thus, we can extract Af~ 
by subtracting the measured asymmetry for Df — > 
from that for D° K~ir + . An A^ correction map 
is obtained as follows; N£^ K ' 7r+ and N°^ K+ *~ are 
corrected according to the reconstructed asymmetry for 
D+ ^7r+ in bins of (p£ h , cos 9^ b , cos6»g MS ). Sub- 
sequently, corrected N°"^ K ~ n+ and N^ K+ ^ val- 
ues are determined in bins of momentum and po- 
lar angle in the laboratory frame, (p*fe, cosO^X)- From 
the corrected values of N£°^ K 



r+ and Np c ° c ^ K+ ^ 



bins of {p l ^f, cos 9 l ^%), we obtain an Af- map that is 
used to correct Af + [1(3] for the KgK + final states. By 
subtracting from the reconstructed asymmetry of 

Df-. — > KgK + , we obtain the corrected reconstruction 



asymmetry Af + for Df, K^K + : 



1 (a 

ircc 



= A T , 

— A {a) 

— FB 



A 



K+ 



D + ,^K°„K+ 
A <s) 
A CP 



(4) 



As shown in Eq. (|4J), A rC c includes not only 

an Acp component but also an App component. Since 
Aqp is independent of all kinematic variables, while App 

is an odd function of cos6>^ s , and thus vanishes when 

(») 

n+^K° s K+ 
integrated over it, we measure j4 re c 



as a func- 
tion of cos rp ne a c p component in Eq. ^ is then 

(«) 

extracted according to Eq. (j5aj) , using the above symme- 
try properties. We also extract the App component in 
Eq. (0} using Eq. (|5b|). 



A (s) 
A CP 



= [^'"^^"(COS^MS 



1 D « 

irec 



(-cos^ s )]/2, 



(5a) 



.4 



FB 



[A 



ICC 



(COS^MS) 



A (s) 



(-COS^MS )]/2 _ 



(5b) 



Figure [3] shows A 



tion of |cos6^ s | 



the |cos^MS| 



CP 



>K° S K+ D+ - 

and App 



as a func- 
Calculating a weighted average over 



bins, we obtain A 



CP 



(-0.16 ± 



0.58)% and A^ KsK = (+0.12 ± 0.36)% where the 
uncertainties are statistical only. The observed App val- 
ues decrease with cos #cms ag ex pected from the leading- 



D 



00 



order prediction [llj. The observed deviations from the 
prediction are expected due to higher order corrections, 
and are in agreement with the measured asymmetries in 
the K + K~ and ir + ir~ final states 

The dominant source of systematic uncertainty in the 

D t)^ K S K+ . K - 

A C p measurement is the uncertainty in Ay , 

which has several sources: the systematic uncertainty 

in the A rC c measurement (0.18%); statistics of the 

D° — > K~ir + sample (0.06%); the systematic uncertainty 
due to the choice of binning for the 2-D map of A*f 
(0.04%); and a possible A CP in the D° -> R-ir+ fi- 
nal state from the interference between decays with and 
without D° — D° mixing. The latter uncertainty is esti- 
mated from the 95% confidence level upper limit on the 
CP violating asymmetry, Aqp = —y sin S sin lp^/R [l3j ]. 
using the world average of D° — D° mixing and CP vio- 
lation parameters and is found to be 0.01%. We also 

binnings to estimate the sys- 



considcr different cos 



gCMS 



tematic uncertainty due to the choice of cos 6>™ s binning 
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TABLE I: Summary of systematic uncertainties. aAcp is the systematic uncertainty in Acp- 



Source 






aA CP s (%) 


oA CP s (%) 


aA CP s (%) 




— y (f>Tv + statistics 


0.18 


0.18 








A Ic i binning 
M(K + K~) window 


0.03 
0.03 


0.03 
0.03 








Dt 4>' K+ statistics 


- 


- 


0.18 


0.18 




A rc i binning 
M(K + K~) window 


- 
- 


- 
- 


0.03 
0.03 


0.03 
0.03 


D° — > K~n + statistics 
A^ binning 
Possible Agl^^ 


- 


- 


0.06 

a a a 
U.04 

0.01 


0.06 

A A/1 

0.U4 
0.01 


cos#^+ B binning 








0.06 


0.06 


Fitting 




0.04 


0.27 


0.12 


0.05 


^"/^-material effects 




0.06 


0.06 


0.06 


0.06 


Total 




0.20 


0.33 


0.25 


0.22 



(0.06%). Systematic uncertainties due to the fitting pro- 
cedure and _ftT°/_ft' -material effects are described above 
and included in Table U where the total systematic un- 
certainties of the Acp measurements are summarized. 
Combining all systematic uncertainties in quadrature, 

we obtain A^ + p ^ K '^ K+ = (-0.16 ± 0.58 ± 0.25)% and 

Acp^ K ° sK+ = (+0.12±0.36±0.22)% where the first un- 
certainties are statistical and the second are systematic. 
Table [TT] summarizes our results, present best measure- 
ments and expected Acp from K° — K° mixing Q. 



TABLE 11: Summary of the Acp measurements. The first 
uncertainties in the second and third columns are statistical 
and the second are systematic. DCS decay contributions are 
ignored for the decays denoted by t's in the fourth column. 

Belle (%) Rcf. [15] (%) Ref. [3] (%) 

Acp^ K ^ + -0.71±0.19±0.20 -1.3±0.7±0.3 -0.332 1 " 
Act^ K ° sn+ +5.45±2.50±0.33 +16.3±7.3±0.3 +0.332 
Acp^ K s K+ -0.16±0.58±0.25 -0.2±1.5±0.9 -0.332 
A^%~^ K ^ K+ +0.12±0.36±0.22 +4.7±1.8±0.9 -0.332+ 
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FIG. 3: Measured Acp and Afb values for D 



« 



K° S K+ as 



a function of | cos #^ms | . The dashed curves show the leading- 



00 



order prediction for A c p B . 

In summary, with a 673 fb _1 data sample collected 
with the Belle detector at the KEKB asymmetric-energy 



we have searched for CP violation in 
the charged charmed meson decays Dt-. — > K®ir + and 



D 



0) 



No evidence for CP violation is ob- 



served. Our results are consistent with the SM (see Ta- 
ble [n| and provide the most stringent constraints to date 
on models for beyond the SM CP violation in these de- 
cays y|. 
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